Vortex state in Bose-Fermi mixture with attraction between bosons and fermions 
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Vortex states in the mixture of ultracold atomic clouds of bosons and fermions are investigated 
using the effective Hamiltonian for the Bose subsystem. A stability of the Bose system in the case 
of attractive interaction between components is studied in the framework of variational Bose wave 
function and Thomas-Fermi approximation. It is shown that the critical number of bosons increases 
in the presence of the vortex. 
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I. INTRODUCTION 

Since the first realization of Bose-Einstein condensa- 
tion (BEC) in ultracold atomic gas clouds Q, H , stud- 
ies in this direction have yielded unprecedented insight 
into the quantum statistical properties of matter. Be- 
sides the studies using the bosonic atoms, growing in- 
terest is focused on the cooling of fermionic atoms to a 
temperature regime where quantum effects dominate the 
properties of the gas yj. 5, 6J. This interest is mainly mo- 
tivated by the quest for the crossover between a BEC and 
Bardeen-Cooper-Schrieffer(BCS) superfluid in ultarcold 
atomic Fermi gasses 0, H, ■ 

Strong s-wave interactions that facilitate evapora- 
tive cooling of bosons are absent among spin-polarized 
fermions due to the exclusion Pauli principle. So the 
fermions are cooled to degeneracy through the mediation 
of fermions in another spin state 0, Ei or via 
a buffer gas of bosons d El El (sympathetic cooling) . 
The Bose gas, which can be cooled evaporatively, is used 
as a coolant, the fermionic system being in thermal equi- 
librium with the cold Bose gas through boson-fermion 
interaction in the region of overlapping of the systems. 

However, the physical properties of Bose-Fermi mix- 
tures are interesting in their own rights and are the sub- 
ject of intensive investigations including the analysis of 
ground state properties, stability, effective Fermi-Fermi 
interaction mediated by the bosons, and new quantum 
phases in an optical lattices El El El El 111 El 
Several successful attempts to trap and cool mixtures of 
bosons and fermions were reported. Quantum degener- 
acy was first reached with mixtures of bosonic 7 Li and 
fermionic 6 Li atoms |g, llCfl. Later, experiments to cool 
mixtures of 23 Na and 6 Li |l9(. as well as 87 Rb and 40 K 
[Tl| . to ultralow temperatures succeeded. 

Although Bose-Einstein condensation and supcfluidity 
are closely connected, they do not occur together in all 
cases. For example, in lower dimensions one can observe 
superfluid without BEC. In principle, rotational proper- 
ties of a Bose-Fermi mixture could directly reveal super- 
fluidity in such systems. Quantized vortices in a rotating 



gas provide direct conclusive evidence for superfluidity 
because they are a consequence of of the existence of a 
macroscopic wave function that describes the superfluid. 
Recently such experimental evidence was obtained for the 
superfluidity in strongly interacting Fermi gas |l2j . 

In this article we study the instability and collapses 
of the trapped boson-fermion mixture due to the boson- 
fcrmion attractive interaction in the presence of the quan- 
tized vortices , u sin g the effective Hamiltonian for the 
Bose system El EH- We analyze quantitatively prop- 
erties of the 87 Rb and 40 K mixture with an attractive 
interaction between bosons and fermions. The stability 
of this system without vortices was recently studied E9 j 
and good agreement with experiment by Modugno and 
co-workers was found. As was shown in the experi- 
ment , as the number of bosons is increased there is 
an instability value Nb c at which a discontinues leakage 
of the bosons and fermions occurs, and collapse of boson 
and fermion clouds is observed. In this article we esti- 
mated the instability boson number Nb c for the collapse 
transition in the presence of the vortices as a function of 
the fermion number and temperature. 



II. EFFECTIVE BOSE HAMILTONIAN 

First of all we briefly discuss the effective boson Hamil- 
tonian [TM Eq . Our starting point is the functional- 
integral representation of the grand-canonical partition 
function of the Bose-Fermi mixture. It has the form 
ElEHH: 



Z 



D[cj>*}D[cj>}D[r]D[^]e^^~{SB{rA)MX) 

+ s F {r^)+s int {4>*^,r^))}- 

and consists of an integration over a complex field </>(t, r) , 
which is periodic on the imaginary-time interval [0, H/3], 
and over the Grassmann field ^(t, r), which is antiperi- 
odic on this interval. Therefore, </>(r, r) describes the 
Bose component of the mixture, whereas ip(r,r) corre- 
sponds to the Fermi component. The term describing 
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the Bose gas has the form: 

eh/3 



where 



d h 2 v 2 



2m B 



S B {<jy\<i>) = J drj dvL*(r,v) (h^- 

+ V B (r)- AlB )0(r ; r) + ^|^(r,r)| 4 }. (2) 

Because the Pauli principle forbids s-wave scattering be- 
tween fermionic atoms in the same hyperfine state, the 
Fermi-gas term can be written in the form: 



+ y F (r)-// J? )V(r,r)}. 



a n 2 v 2 



2mp 



(3) 



The term describing the interaction between the two 
components of the Fcrmi-Bosc mixture is: 

/•ft/3 r 

S int (cf>*,<i>,r,i>)=9BF I drj dr|^(r,r)| 2 |0(r,r)| 2 , 

(4) 

where g B = Ai:h 2 a B /m B and g B p = 2i:h 2 a B p /mi, 
mi = rn B m,p I (m B + mp), m B and mp are the masses 
of bosonic and fermionic atoms respectively, a B and a B p 
are the s wave scattering lengths of boson-boson and 
boson-fermion interactions. 

Integral over Fermi fields is Gassian, we can calcu- 
late this integral and obtain the partition function of the 
Fermi system as a functional of Bose field </>(r, r) . It may 
be shown [TBI HtH that the effective bosonic Hamiltonian 
may be written in the form: 



H, 



eff 



dr(^—\V(b\ 2 + (Vb(t) - mb)I0| 2 + 



[2m B 

+ 9 -fW + feff{\4>\)}. (5) 

where 

feff = ~^- X I V~edeln(l + e^ 
e 3 / 2 de 



) = 



i + ^(e-A) ' 



(6) 



and e = p 2 /2mp, jl = up — Vf{y) — g B p\4>(T,r)\ 2 and 

k = 2 1 / 2 m 3 / 2 /(3^ 2 ft 3 ). 

The first three terms in J5| have the conventional 
Gross-Pitaevskii [24| form, and the last term is a result 
of boson-fermion interaction. In low temperature limit 
p/(k B T) 3> 1 one can write feff{\4>\) m the form: 

feffM) = ~l^ 5/2 - ^n{k B T) 2 ~^ 2 . (7) 
As usual, up can be determined from the equation 



N F = drn F (r). 



(8) 
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At low temperatures we have: 



n F (r) 



ntf' 2 
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7T K 
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(k B T) 2 . 



(9) 



(10) 



In general case the Bose and Fermi systems have dif- 
ferent temperature scales, and Eqs. Q and (|10|) may be 
useful for studying the temperature behavior of Bose sys- 
tem, including the calculation of the critical temperature. 
For example, the characteristic temperature for the Bose 
system - the transition temperature for the ideal Bose gas 
- is 0: k B T° = 0MHuj b (\Nr) 1/3 . The Fermi temper- 
ature for a pure system is [23] k B T F = hu!p(6\N F ) 1 ^ 3 . 
Taking into account that uip = ^/m B jmpU) B , one can see 
that for m B > mp and approximately the same numbers 
of bosons and fermions one can safely use Eqs. J7J) and 
(|10|l to describe the behavior of the Bose system. 

Let us consider now the 87 Rb and 40 K mixture with an 
attractive interaction between bosons and fermions |ll| . 
The parameters of the system are the following: a B = 



5.25 nm, a B p — —21.7 



+4.3 
-4.8 



K and Rb atoms were 



prepared in the doubly polarized states \F — 9/2,mp = 
9/2 > and |2,2 >, respectively. The magnetic poten- 
tial had an elongated symmetry, with harmonic oscilla- 
tion frequencies for Rb atoms uj Bjr = uj b = 2ir x 215-ffz 
and uj B . z = Xlo b = 2tt x l6.3Hz,ujp = \/m B /mpu B rj 
\A1ui b , so that m B Lo 2 B /2 = mpUJp/2 = V~o- The collapse 
was found for the following critical numbers of bosons 
and fermions: N Bc « 10 5 ; N K « 2 X 10 4 . 

At the zero temperature limit, expanding f e ff{\4>\) up 
to the third order in g B p we obtain the effective Hamil- 
tonian in the form: 



H, 



eff 



= / dr 



2m B 
9eff 1 ,,a h- 



■H 4 



+ (V eff (r)-n B )\^ 



9%fW 



(11) 



where 



Veff(r) = ( 1 - ~*A#W) \m B u? B (p 2 + \ 2 z 2 \\2) 



° 1/2 2 /-,Q\ 

9eff = 9B - 2«A*f 5bf> (I 3 ) 

and p 2 = x 2 + y 2 . 

In principle, one can study the properties of a Bose- 
Fermi mixture with the help of f e f / without any ex- 
pansion. However, the form of the Hamiltonian lllfl gives 
the possibility to get a clear insight into the physics of the 
influence of the Fermi system on the Bose one (see discus- 
sion below). It may be easily verified that the expansion 
of the function f(x) = (1 + x) b l 2 (see Eq. Q) up to 
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the third order in x gives a reasonably good approxima- 
tion for f(x) even for rather large values of x, in contrast 
with the higher order expansions, so one can safely use 
Eq. as a starting point for the investigation of the 
properties of the Bose subsystem. 

In derivation of Eqs. (|1 111 3|> we also use the fact that 
due to the Pauli principle (quantum pressure) the radius 
of the Bose condensate is much less than the radius of the 
Fermi cloud Rf ~ y Hf/Vq, so one can use an expansions 
in powers of Vf{y)/ fip. 

From Eq. I|12|) one can see that the interaction with 
Fermi gas leads to modification of the trapping potential. 
For the attractive fermion-boson interaction the system 
should behave as if it was confined in a magnetic trapping 
potential with larger frequencies than the actual ones, in 
agreement with experiment [Tl| . Boson-fermion interac- 
tion also induces the additional attraction between Bose 
atoms which does not depend on the sign of gsF- 

The last term in i? e // (lllfl corresponds to the three- 
particle elastic collisions induced by the boson-fermion 
interaction. In contrast with inelastic 3-body collisions 
which result in the recombination and removing parti- 
cles from the system [3(J, this term for qbf < leads 
to increase of the gas density in the center of the trap in 
order to lower the total energy. The positive zero point 
energy and boson-boson repulsion energy (the first two 
terms in Eq. Hllfl ) stabilize the system. However, if the 
central density grows too much, the kinetic energy and 
boson-boson repulsion are no longer able to prevent the 
collapse of the gas. Likewise the case of Bose condensate 
with attraction (see, for example, 0H3>0|)i the col- 
lapse is expected to occur when the number of particles 
in the condensate exceeds the critical value N Be- 



lli. VARIATIONAL APPROACH 



The qualitative (and even quantitative) picture of the 
stability of the system can be obtained in the framework 
of the variational approach. In the system without vortex 
the critical number Nb c can be calculated using the well- 
known ansatz for the Bosonic wave function 24]: 



(r) 



N B X 



w 3 a 3 ir 3 / 2 



1/2 



exp 



(p 2 + A 2 z 2 )' 



2w 2 a 2 



(14) 



where w is a dimensionless variational parameter which 
fixes the width of the condensate and a — y^h/mBOJB- 

In this case the variational energy Eb has the form: 



Eb 



A 1 



NsfujJB 



bw 2 



ciN B , c 2 N 2 



(15) 



1 - 



V2, 



-Kii F g B F 



C2 



2 \ 9B - ^Vf 9bf 



1/2 9l 
A 2 



A 



(27r) 3 / 2 fiw B o 3 ' 



8fi y 29BF 3 3 / 2 TT 3 fuv B a e ' 



This energy is plotted in Fig.l as a function of w for 
several values of Nb- It is seen that when Nb < Nb c 
there is a local minimum of Eb which correspond to a 
metastable state of the system. This minimum arises due 
the competition between the positive first three terms in 
Eq. (|15|) and negative fourth term. The local minimum 
disappears when the number of bosons Nb exceeds the 
critical value which can be calculated by requiring that 
the first and second derivatives of Eb vanish at the crit- 
ical point. In this case the behavior of Eb is mainly 
determined by the second and fourth terms in Eq. Ijl5|l . 
For Nk — 2 x 10 4 and asF = — 19.44nm we obtain 
Nbc ~ 9 x 10 4 in a good agreement with the experiment 
jllj . It is interesting to note that the critical number of 
Bose atoms in Bose-Fermi mixture is about two orders 
larger than the critical number for the condensate with a 
purely attractive interaction. For example, in the exper- 
iments with trapped 7 Li it was found that the critical 
number of bosons is about 1000. 

If a vortex with angular momentum hi is present in the 
system, the variational wave function may be written in 
the form HI: 



r = 



(wa) 5 7r 3 / 2 



pexp 



p 2 + A 2 z 2 
2uj 2 a 2 



(16) 



In this case the variational energy has the form: 



E b 2 + X 2 + 2l 2 1 „ 2 dN B C 2 N 2 



N B huj B 



Buj z 



Ci = 2 Cl ' 
C 2 = -c 2 . 



From Fig. 1 one can see that the stability of the system 
with a vortex is higher (in this case the critical number 
of bosons Nbc ~ 135 x 10 3 ) due to the centrifugal effect 
(compare the first terms in Eqs. Ijl5|l and Ijl6(l and de- 
creasing the number density of the bosons in the center 
of the trap (fourth terms in Eqs. i|15|) and i|16[l . 

In Eq. (fTK|l we use [i° F = hLUFlQXNp] 1 ^ as the chemical 
potential of the Fermi system [ip . The corrections to /ip 
due to interaction with the Bose system have the form: 

Hf = np[l + mi], where toi = l/2KgBFfi% 1 ^ 2 Ns/Np. 
It may be shown that nil « 0.09 for the values of the 
parameters used in these calculations. The correction to 
the chemical potential for the vortex state is the same. 
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FIG. 1: Variational energy 



N B hu 



as a function of w for 



various numbers of bosons. Upper lines correspond to the 
energy of the system with a vortex. 



FIG. 3: Critical angular vortex velocity as a function of a 
logarithm of the number of bosons Nb ■ 




5,5 6,0 6,5 

I09 10 (N F ) 



7,0 





0,72-1 




0,70- 




0,68- 




0,66- 




0,64- 




0,62- 




0,60- 




0,58- 




0,56- 




0,54- 




4,2 4,4 4,6 4,8 5,0 5,2 5,4 5,6 5,8 6,0 

l°9 10 (N F ) 



FIG. 2: Critical number of bosons Nbc as a function of the 
number of fermions Nf at T = 0. 



FIG. 4: Critical angular vortex velocity as a function of a 
logarithm of the number of fermions Nf- 



Upon increasing the number of fermions, the repulsion 
between bosons decreases leading to the collapse for the 
smaller numbers of the bosonic atoms. In Fig. 2 the crit- 
ical number of bosons Nb c is represented as a function of 
the number of fermions for both cases: the state without 
vortex and the state with the vortex. 

It should be noted that the critical number of bosons 
Nbc is extremely sensitive to the precise value of the 



boson-fermion s wave scattering length. This is illus- 
trated in Fig. 3 in 28]. 

The critical velocity of the vortex state can be found 
from the formula H^: f2 c = (E b {l = 1) - E b )/hN B - In 
Figs. 4 and 5 f2 c is shown as a function of ./Vs and Nf- 

The curves in Figs. 3 and 4 are ending at the critical 
points which correspond to the critical numbers of bosons 
and fermions at which the collapse of the system occurs. 
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The minima in Fig. 4 and intersections of the curves in 
Fig. 5 seem to be artifacts of variational approach to 
calculation of the critical angular velocity in the vicinity 
of the critical points. It should be noted that the num- 
ber of bosons in the system is rather large. In this case 
the qualitatively and even quantitatively correct results 
may be obtained in the framework of the Thomas-Fermi 
approximation (TFA). 



IV. THOMAS-FERMI APPROXIMATION 

In the Thomas- Fermi approximation the kinetic energy 
terms are ignored. It has been shown that in the case of 
one component condensates the TFA results agree well 
with the numerical calculations for large particle num- 
bers, except for a small region near the boundary of the 
condensate 0, 12H |2|j. In fact, even for a small num- 
ber of particles the TFA still usually gives qualitatively 
correct results. The TFA provides an excellent starting 
point of study of the vortex states in Bose condensates 
(see, for example, pi I2I ) . 

The Gross-Pitaevski equation that follows from the ef- 
fective hamiltonian (|llfl : 
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FIG. 5: Critical number of Bosons Nbc as a function of the 
number of Fermions Nf in TFA. 
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A0 + (Krff - Ms) + 5efl#| 2 + 



+ 



3 k 



(17) 



(18) 



In TFA one can neglect the kinetic energy. In this case 
boson density has the form: 



5cff 



I0I 2 = e{(i B - v eS )- 



1 + Jl + (Mi 



y s 3ng% p 



(19) 



In the limit qbf — ► one has the conventional expres- 
sion for the boson density [2~H 12^: 

G"B - Ves) 



\0\ 2 = 0(p B -V eS ) 



9oS 



Equation (|19|) for the boson density n(r) = |</>(r)| 2 may 
be rewritten in the form: 



Here 



n = r t /a h ; R 2 



R 2 



2MB 



corns ui B a\ ' 



'H u / o „ 3 •> ^max 



, --max o 3 2 2' 

3 ng BF 3 ng BF c m B uj B a z h 
3 
2' 



1 1 ^ !/ 2 \ 
CQ = I 1 - ^K^p 9BF ■ 



Rmax is the radius of the condensate. The critical num- 
ber of bosons can be determined from the condition that 
the expression under the square root in l|20|l is positive. 
In Fig. 5 the critical number of bosons is shown as a 
function of the number fermions (lower curve). 

Let us now consider a trap rotating with frequency 
fl along the z-axis. For vortex excitation with angular 
momentum hi , the condensate wave function is given by 



Mr) = VM^je 11 *. 



(21) 



In a rotating frame the energy functional of the system 
is 



E rot (l) = E{cj>i) + / cPrifflitXldyifr) 



(22) 



After substituting the wave function for the vortex ex- 
citation (|21|) in Eq. I|22l) . the effective confinement po- 
tential for the bosons becomes l 2 h 2 /2mp 2 + V, where 
V = muj(p 2 + X 2 z 2 )/2 and p 2 = x 2 + y 2 . So within the 
TFA the density of the vortex state has the form: 



n,(r) =n(0) 1- 



X 2 z 2 



cop 2 



R 2 



R 2 



(23) 



The important new qualitative feature of a vortex in 
the TFA is the appearance of a small hole of radius £, 
£ cx I /R m ax, but the remainder of the condensate density 
is essentially unchanged. The fractional change in the 
chemical potentials caused b y th e vortex (p'(l) — p')/p' 
can be shown to be small |25ll27| . of the order of 1/N 4 / 5 . 
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In the calculation of physical quantities involving the con- 
densate density it is sufficient to retain the no- vortex den- 
sity and simply cut off any divergent radial integrals at 
the appropriate core sizes £ = l/R ma x- Note that us- 
ing the unperturbed density for calculation of the vortex 
properties corresponds to the hydrodynamic limit. 

In Fig. 6 the density distribution of bosons calculated 
from Eqs. (|2(J|I and (|23|l is shown as a function of radius 
for different numbers of bosons and fermions. 

The critical number of bosons in the presence of a vor- 
tex can be calculated from the condition that the ex- 
pression under the square root in (|23l) is positive. This 
number is shown in Fig. 5 (the upper curve) as a function 
of the number of fermions. As was expected the presence 
of the vortex increases the stability of the system against 
the collapse transition. 

From Eqs. (|22l2d|) and (|llli>0|l and Eq. fl c = (E b (l = 
1) — Eb)/hNB one can find the critical angular velocity 
in TEA. The critical angular velocity as a function of the 
numbers of bosons and fermions correspondingly is shown 
in Figs. 7 and 8. One can see that the results are free 
from the mistakes which are present in the variational 
approach (see Figs. 3 and 4). 

V. FINAL REMARKS 

Finally, we make a short remark on the nature of the 
collapse transition. In this article we found the insta- 
bility point of the Bose-Fermi mixture with attractive 
interaction between components with and without vor- 
tex. A strong rise of density of bosons and fermions (see 
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0,48- 

0,46 J , — 

100000 



FIG. 8: Critical angular vortex velocity as a function of a 
logarithm of the number of fermions Nf in TFA. 

Eq. HI U| l) m the collapsing condensate enhances intrin- 
sic inelastic processes, in particular, the recombination in 
3-body interatomic collisions, as is the case for the well- 
known 7 Li condensates |30j . In the presence of a vor- 
tex there appears a hole in the middle of the condensate. 
This reduces the maximum density of the condensate (see 
Fig. 6) and increases the critical number of bosons. Re- 
cently M. Yu. Kagan and coworkers suggested the new 
microscopic mechanism of removing atoms from the sys- 
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tern which is specific for the Bose-Fermi mixtures with 
attraction between components and is based on the for- 
mation of the boson- fermion bound states [^2] • It seems 
that the description of the evolution of the collapsing 
condensate should include both these mechanisms. 



VI. ACKNOWLEDGEMENT 

VNR acknowledges valuable discussions with M. Yu. 
Kagan. The work was supported in part by the Russian 



Foundation for Basic Research (Grant No 05-02-17280) 
and NWO-RFBR Grant No 047.016.001. STC is partly 
supported by a grant from NASA. NMC acknowledges 
the support from RFBR (Grant No 04-02-08159), Gov- 
ernment Grant -4611.2004.2 and Russian Science Sup- 
port Foundation. 



[16] 



M. N. Anderson, J. R. Ensher, M. R. Matthews et al., Sci- [17 
ence 269, 198 (1995). 
K.B. Davis, M.-O. Mewes, M. R. Andrews et al, Phys. [18 
Rev. Lett. 75, 3969 (1995). 

C. C. Bradley, C. A. Sackett, and R. G. Hulet, Phys. Rev. [19 
Lett. 78, 985 (1997). 

B. DeMarco and D. S. Jin, Science 285, 1703 (1999). [20 
A. G. Truscott, K. E. Strecker, W. I. McAlexander et al, 
Science 291, 2570 (2001). [21 
K.M. O'Hara, S.L. Hemmer, M. E. Gehm et al, Science 

298, 2179 (2002). [22 

C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Nature [23 
424, 47 (2003). 

M. Greiner, C. A. Regal, and D. S. Jin, Nature 426, 537 [24 

(2003) . 

S. Jochim, M. Bartenstein, A. Altmeyer et al, Science [25 
302, 2101 (2003). 

F. Schreck, L. Khaykovich, K. L. Corwin et al, Phys. Rev. [26 
Lett. 87, 080403 (2001). [27 

G. Modungo, G. Roati, F. Riboli et al, Science 297, 2240 [28 
(2002). 

M. W. Zwierlein, J. R. Abo-Shaeer, A. Schirotzek, [29 
C.H. Schunck, and W. Ketterle, Nature, 435, 1047 
(2005). [30 
T. Miyakawa, T. Suzuki, and H. Yabu, Phys. Rev. A 64, 
033611 (2001). [31 
R. Roth, Phys. Rev. A 66, 013614 (2002). 

S.T. Chui and V.N. Ryzhov, Phys. Rev. A 69, 043607 [32 

(2004) . 

S.T. Chui, V.N. Ryzhov, and E. E. Tareyeva, JETP 
Lettes 80, 274 (2004). 



M. J. Bijlsma, B. A. Heringa, and H. T. C. Stoof, Phys. 
Rev. A 61, 053601 (2000). 

M. Lewenstein, L. Santos, M. A. Baranov, and H. 
Fehrmann, Phys. Rev. Lett. 92, 050401 (2004). 
Z. Hadzibabic, C. A. Stan, K. Dieckmann et al, Phys. 
Rev. Lett. 88, 160401 (2002). 

J. Goldwin, S. B. Papp, B. DeMarco, and D. S. Jin, Phys. 
Rev. A 65, 021402(R) (2002). 

V. N. Popov, Functional Integrals in Quantum Field The- 
ory and Statistical Physics (Reidel, Dordrecht, 1983). 
H.T. C. Stoof, cond-matt/9910441. 

D.A. Butts and D.S. Rokhsar, Phys. Rev. A 55, 4346 
(1997). 

F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999). 

A. L. Fetter and A. A. Svidzinsky, J. Phys.: Condensed 
Matter 13, R135 (2001). 

G. Baym and C.J. Pethick, Phys. Rev. Lett. 76, 6 (1996). 
S. Sinha, Phys. Rev. A 55, 4325 (1997). 

S. T. Chui, V. N. Ryzhov and E. E. Tareyeva, JETP 
91,No 6, pp. 1183-1189 (2000). 

S. T. Chui, V. N. Ryzhov, and E. E. Tareyeva, Phys. 
Rev. A 63, 023605 (2001). 

Yu. Kagan, A. E. Muryshev, and G. V. Shlyapnikov, 
Phys. Rev. Lett. 81, 933 (1998). 

Yu. Kagan, G. V. Shlyapnikov, and J. T. W. Walraven, 
Phys. Rev. Lett. 76, 2670 (1996). 

M.Yu. Kagan, I. V. Brodsky, D.V. Efremov, A.V. 
Klaptsov, cond-mat/0209481. 



